Band structure model of magnetic coupling in semiconductors 
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We present a unified band structure model to explain magnetic ordering in Mn-doped semicon- 
ductors. This model is based on the p-d and d-d level repulsions between the Mn ions and host 
elements and can successfully explain magnetic ordering observed in all Mn doped IF VI and III- 
V semiconductors such as CdTe, GaAs, ZnO, and GaN. This model, therefore, provides a simple 
guideline for future band structure engineering of magnetic semiconductors. 
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Mn-doped II- VI and III-V diluted magnetic semi- 
conductors (DMS) have many unique magneto-optical, 
magneto-electrical, and magneto-transport properties 
that are essential for future-generation spintronic device 
applications (H Q- [J Dl These materials also present 
many interesting behaviors. For example, Mn-doped 
II- VI semiconductors generally have antiferromagnetic 
(AFM) ground state, whereas Mn doped III-V semicon- 
ductors mostly have ferromagnetic (FM) ground state. 
Of great interest, some of the systems, such as Mn-doped 
GaN, can be either FM or AFM, depending on the Mn 
concentration, carrier density, or pressure It will 
be, therefore, rewarding to understand the mechanism 
of magnetic coupling in these systems in order to design 
functional spintronic materials. 

The mechanism that is responsible for the magnetic 
coupling in Mn doped III-V semiconductors has been 
widely discussed over the last few years 0, ED ■ Sev- 
eral models have been proposed to explain the phenom- 
ena, including the Zener/RKKY and, super- and double- 
exchange models. Although these models are quite suc- 
cessful in explaining some of the systems, they often lack 
universality and transparency, and are difficult to com- 
pare directly with ab initio band structure calculations. 
In this paper, using a band structure approach and level 
repulsion model, we will describe a unified picture to 
explain magnetic ordering in Mn-doped III-V and II- VI 
semiconductors and relate them to the previous models. 
We show that our model can successfully explain mag- 
netic ordering observed in all Mn-doped II- VI and III- 
V semiconductors such as CdTe, GaAs, ZnO, and GaN. 
Therefore, it provides simple guidelines for future band 
structure engineering of magnetic semiconductors. 

Our model is derived from the coupling and level repul- 
sion between the magnetic ions d and the host element's p 
states. When Mn atoms are incorporated into a semicon- 
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FIG. 1: Schematic model for p-d and d-d couplings, (a) p-d 
coupling with ti v above tid- (b) p-d coupling with t2 P below 
t2d- (c) d-d coupling between the tid states with majority spin 
level coupled to majority spin level and minority spin level 
coupled to minority spin level, and (d) d-d coupling between 
the t2d states with majority spin level coupled to minority 
spin level, (e) Same as (c), but between the ed state, (f) 
Same as (d), but between the ed state. 



ductor, they will introduce d levels either inside the va- 
lence band or above the valence band maximum (VBM) 
of the semiconductor, depending on the host material. 
For example, for Mn in GaAs and CdTe, the occupied 
majority d bands are below the host VBM, which con- 
tains mostly anion p states. The unoccupied minority d 
band is above the VBM. On the other hand, in ZnO and 
GaN, the d levels are above the VBM. The introduced 
d orbitals will interact with the host p states forming 
hybrid pd orbitals. The hybrid orbitals will then also in- 
teract with each other. The hybridization occurs because 
the Mn d states in a thetrahedral crystal field split into 
t-2 and e states, whereas the VBM also has the £2 symmc- 
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try, with large anion p-character. The two states with the 
same t 2 symmetry can couple strongly with each other. 
Furthermore, without spin-orbit (SO) coupling, only the 
states with the same spin configuration can couple to 
each other 0, 0] . Figure 1 shows the possible scenar- 
ios for the states that may couple to each other. In Fig. 
la, the host t 2p state is above the magnetic ion t 2 d state, 
and the coupling leads to a level repulsion that pushes 
up the t 2p state by A pd and pushes down the t 2 d state 
by A pd . In Fig. lb, the t 2 d state is above the t 2p state, 
and they are pushed up and down respectively by A pd . 
Figure lc shows the couplings between the two majority 
spin d levels and the two minority spin d levels, which 
lead to splittings 2A dd and 2A dd , respectively. Figure 
Id shows the coupling between the majority spin d state 
and the minority spin d state, which pushes up the mi- 
nority t 2 d state by A^ 2 and pushes down the majority 
spin state by the same amount. The coupling between 
the ed states is similar to that of between the t 2 d states, 
and they are described in Figs, le and If. In principle, 
the strenght of the interactions in Fig. 1 depends on the 
distance and orientation of the Mn pair. From these fig- 
ures, we can see that there will be no energy gain if the 
two coupled states are fully occupied or fully empty. A 
magnetic state can be stabilized only if there are both 
holes and electrons in the coupled states. 

In general, the p-d coupling is larger than the d-d cou- 
pling because the host p orbital is more delocalized and 
Mn d and anion p are nearest neighbors. Within the d- 
d coupling, the coupling between the ed state is much 
weaker than the coupling between the t 2 d state, because 
the ed state is very localized, whereas the t 2 d state is 
less localized due to the coupling with the host p states. 
Because the d-d coupling between the e<j state is qualita- 
tively similar to that of between the t 2 d state, they are 
not discussed explicitly in our model. Furthermore, A dd 
is smaller than A dd when the exchange splitting Cdd is 
not zero because the level repulsion is larger if the en- 
ergy separation of the two coupling states is small. 

With the understanding above, in the following, we will 
discuss the mechanism of magnetic ordering observed in 
Mn doped II- VI and III-V semiconductors. In the discus- 
sion, we will use the two-band coupling model described 
in Fig. 1. First we will discuss the interaction between 
the p and d levels, and then the interaction between the 
d levels. We will assume the coupling is additive. The 
small effect of coupled p-d and d-d interactions will be 
discussed briefly for simplicity, but in general, the treat- 
ment would not change the results qualitatively. 

(i) We discuss first the scenario when the VBM is be- 
tween the Mn majority and minority spin d states. This 
is the case for most Mn-substituted II- VI and III-V sys- 
tems such as CdMnTe and GaMnAs. A schematic plot is 
shown in Fig. 2, where the level repulsions are arranged 
in spin-up and spin-down channels in FM and AFM con- 







Spin -Up 








Spin -Down 




(a) 


(b) 


(c) 


<b') 


(a') 


(a) 


(b) (c) (b') 


(a') 


■j 


,-tto. 








t 2d ' 


,, 


' V 












Ill ■-■■■"::.iaa-' ::::! " 111 
ttt TfT ttt 




1 t,T***-- 




4*1- 




MLt 






b 


•-***- 


-***- 










u, 




4H- 












gnetic 








t M 




, 






.-**».. 










-tr* 




















Anti: 




4*±- 













FIG. 2: Schematic model showing the position of the p and 
d levels and level repulsion between them in FM and AFM 
configurations. In this case, the VBM is in between the ma- 
jority and minority d levels. Note that in (b), (b'), and (c), 
the states have mixed pd characters. 

figurations. From (a) to (b) or (a') to (b'), only the p-d 
coupling is turned on. From (b) and (b') to (c), the d-d 
coupling is turned on. We first look at the FM configu- 
ration. In the spin-up channel, the p-d coupling pushes 
the t 2p state upward by 2A* d , whereas it pushes down 
the t 2d state by the same amount [l^. The net energy 
gain in this process is — 2mi l A 1 i , where rrih is the number 
of holes. The d-d coupling further splits the d levels by 
2A^ d . Because all the t 2 d majority spin levels are fully 
occupied, there is no energy gain in this process. (The p-p 
coupling is already included in the band structure calcu- 
lation for the host, and its effect is the same for FM and 
AFM cases; therefore, it is not discussed here.) In the 
spin-down channel, the p-d coupling lowers the energy of 
the occupied t 2p spin-down state by 2A pd . Because there 
are six electrons in the two t 2p states, the net energy gain 
in this case is — 12A pd . Therefore, the net energy gain for 
the FM configuration is -2m h A pd - \2A 2 pd . In the AFM 
configuration, the situation is the same in the spin-up 
and spin-down channel. When the p-d coupling is turned 
on, the t 2p state is pushed up by the occupied majority 
t 2 d state by A pd , but pushed down by the unoccupied mi- 
nority t 2 d state by A pd . The occupied majority t 2 d state 
is pushed down by A pdl and it is further pushed down 

by the unoccupied minority t 2 d state by A^ 2 . There- 
fore, the net energy gain in the AFM configuration is 
-m h A pd - (12 - m h )A 2 pd - 6A d ' d . Hence, taking into 
account the localization of the host p orbitals around the 
Mn atoms 0], the energy difference between the FM 
and AFM phase is 

AEfm-afm = -am h (Al d + A 2 pd ) + 6A^ 2 , (1) 

where a < 1 decreases when the hole states become more 
localized, and when Mn-Mn distance increases. This re- 
sult suggests that (a) the AFM phase is stabilized by 
the energy 6A^ 2 from coupling between the majority 
and minority spin d states (often denoted as superex- 
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FIG. 3: Schematic model showing the position of the p and 
d levels and level repulsion between them in FM and AFM 
configurations. In this case, the Mn d levels are above the 
VBM. Note that in (b), (b'), and (c), the states have mixed 
pd characters. 



change) [3]. (b) The FM phase is stabilized with energy 
— mh(Ap d + Ap d ) which is proportional to the number 
of holes and the p-d exchange splitting (A pd + A pd ), as 
described in the Zener model. Therefore, to enhance FM 
coupling, one should increase the hole carrier density and 
increase the p-d exchange splitting. However, large p-d 
coupling also leads to large localization of the hole stae, 
thus, a balance between a and A p d is needed; (c) For 
systems where Mn substitution on the cation site does 
not introduce holes (e.g., CdMnTe), the system is always 
more stable in the AFM phase, (d) For systems where 
Mn substitution on the cation site introduces holes (e.g., 
GaMnAs, where each Mn on Ga site introduces one hole), 
because the p-d coupling is larger than the d-d coupling, 
the system in general will have a FM ground state if 
enough holes are present in the system. However, when 
holes are compensated by donor defects, the system can 
revert to the AFM ground state. 

(ii) In the second scenario we discuss the case when 
the Mn d state is above the VBM of the host. This is 
the the case for Mn in ZnO or in GaN. A schematic plot 
is shown in Fig. 3. The analysis is similar to that in 
the first scenario. In the FM configuration and spin-up 
channel, the system not only gains energy through the p-d 
coupling by —2mh,A pd , but also through the d-d coupling 
[lif by —rrihA dd , which put holes at a high energy level 
and electrons at a low energy level. Here, < 3 is the 
number of holes at the tid level, and for simplicity, we 
still use A pd to describe the coupling between the VBM 
and the majority tid state. The energy gain in the spin- 
down channel is — Y2A 2 pd , thus the total energy gain for 
the FM configuration is —m h A dd — 2m h A pd - l2A 2 pd . 
On the other hand, the net energy gain for the AFM 
configuration is -(6 - m h )A dd - 2m h A pd - l2A 2 pd , thus 
the energy difference between the FM and AFM phase is 
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FIG. 4: Total (black) and Mn d projected (red) density of 
states for CdMnTe and GaMnAs: (a) ferromagnetic CdMnTe; 
(b) antiferromagnetic CdMnTe; (c) ferromagnetic GaMnAs; 
and (d) antiferromagnetic GaMnAs. 



These results indicate that when the system has holes 
at the £2d-derived level instead of the pd orbital, (a) the 
stabilization of the FM or AFM phase is not directly re- 
lated to the p-d exchange splitting, but is determined by 
the d — d coupling [J^j terms A dd (often denoted as dou- 
ble exchange) and A dd terms, (b) When the number of 
holes rrih decreases, the AFM phase is stabilized against 
the FM phase. When mu = (e.g., ZnMnO), the system 
is always more stable in the AFM phase, (c) Whether the 
system has an FM or AFM ground state depends not only 
on hole carrier density, but also on the relative strength 
of FM stabilization energy A dd , and the AFM stabiliza- 
tion energy A^ 2 . A^ 2 increases when the exchange split- 
ting edd decreases. This can be achieved if the charge is 
transfered from the majority to minority spin state (e.g., 
from majority tid to minority state), thus reducing the 
magnetic moment and exchange splitting. The majority 
tid state can be pushed upward through increased p-d 
coupling, which can be realized by increasing Mn con- 
centration or applying pressure. Thus, we predict that 
for Mn-doped GaN, when the Mn concentration is high 
or under pressure, the system will have an AFM ground 
state 0. 

To test our models, we performed first-principles total- 
energy calculations for CdMnTe, GaMnAs, ZnMnO, and 
GaMnN. The calculations were performed using an ab 
initio plane wave basis code , based on the local spin 
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FIG. 5: Total (black) and Mn d projected (red) density of 
states for ZnMnO and GaMnN: (a) ferromagnetic ZnMnO; (b) 
antiferromagnetic ZnMnO; (c) ferromagnetic GaMnN; and 
(d) antiferromagnetic GaMnN. 



density functional theory and using ultrasoft pseudopo- 
tentials 0|. For the exchange and correlation poten- 
tial, we used the generalized gradient approximation of 
Perdew and Wang 17]. The Brillouin zone integration 
is performed using the Monkhost-Pack special k points 
scheme 0] of 4x4x4 for the energy differences and 6x6x6 
for the DOS, in a 64 atom supercell with two Mn atoms 
as first fee neighbors. Interactions of further neighbors 
were considered elsewhere 01- We considered the zinc- 
blende alloy and assumed that the same results also hold 
for the alloy in a wurtzite structure. 

Figure 4 shows the total (black) and Mn d projected 
(red) density of states for CdMnTe and GaMnAs in the 
FM and AFM configurations. We see that these two 
systems correspond to scenario (i) where the VBM is be- 
tween the majority and minority Mn d states. For CdM- 
nTe, the substitution of Mn for Cd does not introduce 
holes (Figs. 4a and 4b), therefore, according to our model 
[Eq. (1)], the system should be AFM. For GaMnAs, holes 
are present in the VBM-derived states (Figs. 4c and 4d), 
therefore, our model predicts that it should be FM. In- 
deed, our directly calculated total energy differences be- 
tween FM and AFM configurations for these two systems 
(Table I) agree with those expected from our model. Fig- 
ure 5 shows the total (black) and Mn d projected (red) 
density of states for ZnMnO and GaMnN in the FM and 
AFM configurations. We see that these two systems cor- 



TABLE I: Energy difference between FM and AFM configu- 
rations for Mn-doped semiconductors [x = 6.25%). 
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respond to scenario (ii) where the VBM is below the Mn 
d states. For ZnMnO, no holes are present when Zn is 
replaced by Mn (Figs. 5a and 5b), therefore, according 
to our model [Eq. (2)], the system should be AFM. For 
GaMnN with small Mn conccentration (6.25%), holes are 
created at the Mn 3d-derived level, and the Mn has a 
high-spin configuration (Figs. 5c and 5d), thus our model 
predicts that it should have an FM ground state. Again, 
our direct calculations shown in Table I agree with what 
is expected from our model. 

In summary, we have proposed a model that can suc- 
cessfully explain the stabilization of the ferromagnetic 
or antiferromagnetic phases in Mn-doped II- VI and III- 
V semiconductors. This simple model is based on the 
p-d and d-d level repulsions between the Mn ions and 
host states and can be directly related to band parame- 
ters. Therefore, it should be very useful in understanding 
and engineering diluted magnetic semiconductors with 
desired properties. A similar model could also be devel- 
oped to study other transition metals in semiconductors, 
as long as the position of the d levels and the number of 
holes induced by it were known. 
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